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A theoretical study of the surface energy- loss function of freestanding Pb(lll) thin films is pre- 
sented, starting from the single monolayer case. The calculations are carried applying the linear 
response theory, with inclusion of the electron band structure by means of a first-principles pseu- 
dopotential approach using a supercell scheme. Quantum-size effects on the plasmon modes of the 
thinnest films are found in qualitative agreement with previous work based on the jellium model. 
For thicker films, results show a dispersionless mode at all thicknesses, in agreeement with electron 
energy-loss measurements. For sizeable values of the momentum, the raising of the surface plasmon 
with increasing thickness is retrieved. 
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I. INTRODUCTION 

In thin metallic films, the confinement in the direction 
perpendicular to the film plane gives rise to the quanti- 
zation of the electronic wave functions. As a result of the 
appearance of the so-called quantum-well states (QWS),^ 
the properties of the metallic slabs might strongly depend 
on the exact thickness of the film. This dependence is a 
purely quantum phenomena known as quantum-size ef- 
fect (QSE), which often appears as an oscillatory depen- 
dence of several physical properties on the film thickness. 

Thin lead films exhibit important quantum-size oscilla- 
tions in the layer- by-layer growth,^ first observed by He- 
atom scattering and attributed to interference with the 
quantum- well states. The latter modulate the electron 
density of states (DOS) at the Fermi level (Ep), causing 
oscillations with varying thickness in the superconduct- 
ing critical temperature and the upper critical field, ^ ^ 
interlayer distances,^ island height distributions,'^ zone- 
center phonon frequencies,^^ electronic transport ,'^pho- 
toemission pr opert ies,*^ work functional^' and quasipar- 
ticle lifetimes. ^^ ^^ Also, recently superconductivity was 
discovered in a single lead monolayer (ML) on silicon. ^^ 
Thus, lead films have become an important model sys- 
tem for exploring electronic and structural properties of 
metals on the nanoscale.^^ 

However, to the best of our knowledge, there are few 
experimental studies on the surface dielectric response of 
Pb thin slabs,^^^^^ and no theoretical works. Thus, the 
aim of the present work is performing a computational 
systematic study of the surface energy-loss function of 
Pb(lll) films with varying thickness, starting from the 
single monolayer case, up to the 15 ML thick slab. 

An approximate description of thin film plasmons is 



given by the solution of the Maxwell equations for the 

proper geometry. ^^ It leads to the coupling between the 

classical surface plasmons of the two different surfaces of 

the film. The resulting coupled modes of the film disperse 
as^^«2ij 



UJ± 



^<-'- 



qL\l/2^ 



(1) 



where Up is the bulk plasmon frequency, which is given 
by ujp = ^/4:7^n/m* with n being the average electron 
density and m* the electron effective mass, which in 
terms of the density parameter r^ standing for the av- 
erage inter-electron distance reads as uOp = y3/r|m*. 
The energy splitting between the modes depends on the 
film thickness L and the in-plane momentum transfer q. 
The low-energy mode uj- corresponds to a symmetric in- 
duced charge profile in the direction perpendicular to the 
film plane, whereas the high-energy mode cj+ corresponds 
to an antisymmetric one.^^ As L increases, the coupling 
between the two modes decreases. In the limit L^l/q 
the two film modes are decoupled and the two classical 
surface plasmons of frequency ujp/V2 are retrieved. This 
model ignores the electronic structure of the film. This 
is a serious drawback since the ground state electronic 
structure has been shown to strongly affect the surface 
response to external perturbations. More detailed classi- 
cal models showed the dependence of the surface plasmon 
dispersion on the microscopic details of the surface elec- 
tronic density profile. ^^ ^^ 

On a more quantitative level, the jellium model^^ has 
been used to study the quantum-mechanical electrody- 
namical response of metal slabs,*^^"^ gaining basic in- 
sight into the nature of electronic excitations of metallic 
films. As an example. Yuan and Gao have shown,*^ using 
the jellium model with the electron density correspond- 



ing to Ag, the disappearance of the antisymmetric mode 
ujj^ for g ^ when the film thickness is comparable to 
the Fermi wavelength. Instead, a few discrete interband 
peaks were found. ^^ 

A more precise description of the electron band struc- 
ture in the direction perpendicular to the film planep^'^ 
allowing to describe the surface states which are missing 
in a jellium model, was recently used to study new col- 
lective electronic excitations at metal surfaceP^'^^ and 
thin metal films.'^ Howe ver, t he recipe of the improved 
one-dimensional potentia r^ * ^^ * can not give a satisfactory 
description of the electronic structure of Pb(lll) films. 
Thus, in the present work a first-principles approach is 
used to study the dielectric response of Pb(lll) films. 
Indeed, using an ah initio calculation scheme possible 
anisotropy effects can be studied, which are missing in 
jellium models or in using the potentials of Refs.l28land 
[29^ as they assume in-plane free-electron-like behavior. 

In general, the plasmon modes of the Pb(lll) films are 
found to follow qualitatively the classical dispersion re- 
lation Eq. (fTl). However, for the thinnest slabs QSE are 
found. A dispersionless mode is found at all thicknesses, 
replacing in practice the surface plasmon as the short 
wavelength limit of the symmetric plasmon mode. The 
raising of the surface plasmon with increasing thickness 
is found at short wavelengths. Also, comparison of the 
surface energy-loss function with the momentum trans- 
fer q along different high-symmetry directions showed no 
sizeable anisotropy effects. 

The rest of the paper is organized as follows. In Sec. II 
the details of the ah initio calculation of the surface loss 
function using a supercell scheme are shown. In Sec. Ill 
the calculated ground state electronic structure proper- 
ties are presented, while in Sec. IV the results on the sur- 
face loss function are analyzed in detail. Finally, the main 
conclusions are drawn in Sec. V. Unless otherwise stated, 
atomic units are used throughout, i.e., e^ = ^ = mg = 1. 



II. CALCULATION METHOD 

When a perturbing electric charge is located far from 
one side of the film the differential cross section for its 
scattering with energy uj and in-plane momentum trans- 
fer q is proportional to the imaginary part of the surface 
response function g(q, cj) defined a&^ 



which depends on the film electronic properties only 
{q. = |q|)- This quantity is relevant in the descrip- 
tion of surface collective excitations measured in elec- 
tron energy-loss experiments.' ^^ ' ^^ ' Here x{^ •> ^' •> ^.i ^) is 
the density response function of an interacting electron 
system that determines, within linear response theory, 
the electron density n^"^(2:, q, a;) induced in the system 
by an external potential U®^*(2;, q, cj) according to 



^ind 



(z, q, o;) = I dz'x{z, z^ q, uj)V-^\z' , q, uj). (3) 



g(q,^) = j dz j dz'x{z,z',q^, 



uj)e' 



q{z+z') 



(2) 



The collective electronic excitations in thin films then can 
be traced to the peaks in the surface loss function defined 
as the imaginary part of g, Im[g(q, cj)]. 

In t he fra mework of time-dependent density functional 
theory,'^^'^ X is the solution of the integral equation 



X{z,z',ci,uj) = x°(^,^',q,^)+ / dzi j dz2X^ {z , z , q, uj) 
X [vc{zi,Z2,q)^Kxc{zi,Z2,q,uj)]x{z2,z\q,uj), (4) 



with x^ being the response function of the noninter- 
acting Kohn-Sham electrons. In Eq. ([3| Vc{z^z\q) = 
— — e^l^~^ I stands for the two-dimensional (2D) Fourier 
transform of the bare Coulomb potential and Kxc ac- 
counts for the exchange-correlation (XC) effects. In the 
present work, we use the random-phase approximation 
(RPA) where Kxc is set to zero, i.e., the dynamical 
short-range exchange-correlation effects are ignored. Pre- 
vious studies of collective excitations at the surfacea^^lESl 
and in the bullP^J^ of many metallic systems suggest 
that XC effects should have little impact on the study of 
Pb films. 

For a periodic system, the polarizability can be ex- 
pressed as a matrix in the basis of the reciprocal space 
vectors {G}. As a consequence, Eq. Q becomes a ma- 
trix equation. Then, once the ground state has been ob- 
tained, the starting point of the calculation of the surface 
response function is the evaluation of the matrix elements 
of the noninteracting polarizability 
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where n {rC) is an occupied (unoccupied) band index, k is in the two-dimensional surface Brillouin zone (SBZ), 
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FIG. 1: (Color online) Calculated band structure and density of states of Pb(lll) freestanding films 1 - 3, 6, 10 and 15 ML 
thick, with (solid lines) and without (dashed lines) spin-orbit (SO) coupling in the Hamiltonian. The horizontal dashed-dotted 
lines represent the Fermi level of each film. DOS in arbitrary units. 



/k,n are Fermi factors and £^k,n (^k,n) are Kohn-Sham 
energies (wave functions). Actually in order to speed 
up the calculations, following Refs. "43^ and '45 , first the 
spectral function is calculated and from its knowledge the 
imaginary and real parts of Xq g' ^^^ obtained. Finally, 
the expression for the surface response function in the 
case of a periodically repeated slab reads 

(6) 
Even though only the G = G' = matrix element of 
Xg,G' enters Eq. (|6|, the full three-dimensional (3D) na- 
ture of the polarizability is implicity taken into account 
via the evaluation of Eq. Q as a matrix equation. 

In order to save computational time, Xg g' ^^^ been 
calculated retaining only G = (0,0,^^) reciprocal space 
vectors. Physically, this means that lateral crystal lo- 
cal field effect &^^ were neglected. This approach was al- 
ready found to give indistinguishable results compared 
with the calculations carried using the 3D G's for metal 
surfaces.^ All important 3D effects are included in the 



evaluation of Xg g' through the use of the fully 3D Bloch 
functions and their respective one-electron energies. 

In the present work Pb(lll) films are represented by 
freestanding slabs infinite in the xy plane and periodically 
repeated in the z direction, separated by a vacuum region 
whose thickness here is fixed in all cases as 10 interlayer 
distances of the lead atoms of the film in the z direction. 
Films are not relaxed, representing ideal cuts of the face- 
centered cubic bulk Pb in the (111) direction with the 
bulk experimental lattice parameter of 4.95 A. Thus, the 
in-plane lattice parameter is a = 3.50 A, the interlayer 
distance is c = 2.86 A and the vacuum region thickness 
is d = 28.6 A. However, 4-6 ML thick films were also 
allowed to relax in the z direction and their band struc- 
ture showed small changes compared with their unrelaxed 
counterparts. 

For the density functional theory (DFT) ground state 
calculations, the electron- ion interaction is represented 
by norm-conserving non-local pseudopotentials,'^ and 
the LDA approximation is chosen for the exchange 
and correlation potential, with the use of the Perdew- 



Zungei^^ parametrization of the XC energy of Ceper- 
ley and Alder. ^^ Well- converged results have been found 
with a kinetic energy cut-off of ~220 eV, including from 
- 2200 (1 ML) to - 5300 (15 ML) plane-waves in the 
expansion of the Bloch states. 

For 1-4 ML thick films, the Hamiltonian was also 
solved including the SO term fully self consistently. As a 
centrosymmetric supercell was used in the calculations, 
due to the Kramers degenerac}/^^ the electron energy 
bands are doubly degenerate also when the SO coupling 
is included in the Hamiltonian (see Fig. IT]). 

The calculation of Xq g' (^' ^) ^^^ carried out using a 
Monkhorst Pack 192x192x1 (96x96x1) grid of k vectors 
as the hexagonal SBZ sampling with 3169 (817) k vectors 
in the irreducible part of the SBZ for the 1-5 ML (6-15 
ML) thick films. Up to 500 bands were included in the 
evaluation of Xg g'(^'^) ^^^ ^^^ thicknesses. The width 
of the Gaussian replacing the delta function in the eval- 
uation of Xg g'(^'^) ^^^ ^^^ ^^ ^'^^ ^^' ^ value which 
gave smooth results while not hiding any feature on the 
surface loss function of the films. Well converged results 
are found including 750 plane waves in the expansion of 
the wave functions in the calculation of Xg g' (^' ^) ^^^ 
expanding the size of the polarizability matrix up to 60 
G vectors. 



III. ELECTRONIC STRUCTURE RESULTS 

The results of the ground state calculations showed bi- 
layer oscillations as a function of the slab thickness on 
the density of states at the Fermi level and on the work 
function, with a beating pattern of period 9 ML superim- 
posed (not shown). This is in agreement with previous 
experimental and theoretical studies (see i.e. Refs. ITTI 
and [51]). 

In Fig. [l] the calculated electronic band structure of 
Pb(lll) films of several thicknesses is shown. For a N 
ML thick film, each electron state energy level is un- 
folded in N subbands. The subbands below -6 eV are 
of s character. They are separated by a gap from the 3N 
subbands of p character which form the Fermi surfaces 
of the slabs. As can be seen from Fig. [l] the width of the 
gap is already fixed as ~2 eV for the 3 ML thick film. 

Around the SBZ center (F point) bands show a 
parabolic free-electron-like dispersion. The p bands 
around F present a Pz character, while acquiring an in- 
creasing px^y component as they loss their parabolic-like 
dispersion moving away from F.^^ The Pz states represent 
the QWS of the Pb(lll) films. The present work found 
that the inverse of the energy separation of the QWS 
around Ep is linearly proportional to the film thickness, 
in good agreement with a previous study. ^^ 

As lead is a heavy element (atomic number 82), SO 
interaction has sizeable effects on its energy spectrum. 
As an example, in bulk Pb the SO-induced splitting at 
the BZ center is ~3 eV, and several degeneracies are 
lifted throughout the BZ.I^ In Fig. [l] the band struc- 



ture and DOS for the 1-3 ML thick films is shown 
with (dashed lines) and without (solid lines) SO cou- 
pling included in the Hamiltonian. As can be readily 
seen, SO effects are remarkable only for the single mono- 
layer case, which becomes semimetallic when the SO in- 
teraction is switched on, as a result of the avoiding of 
the band-crossings present for the scalar-relativistic sys- 
tem around the Fermi level. As the slab thickness is 
increased, also the filling of the phase space by the un- 
folding of the subbands increases. Because of the fast 
filling of the phase space, SO effects on the ground state 
of Pb(lll) films become small for slabs as thin as 3 ML 
(see Fig. [l]), as avoiding of the band-crossings is the only 
sizeable SO effect on their band structure. As a conse- 
quence, SO effects are not expected to affect qualitatively 
the films surface loss function, except for the somewhat 
artificial semimetallic single Pb(lll) monolayer. Thus, 
in the present work only scalar-relativistic calculations 
are reported. Note however that a recent first-principles 
study has shown the inclusion of the SO coupling as nec- 
essary in the calculation of the electron-phonon coupling 
and the superconducting temperature of Pb(lll) films.^ 



IV. SURFACE LOSS FUNCTION RESULTS 

A. Isotropy 

In Fig. [2] the calculated surface loss function for the 
3MLs thick Pb(lll) film, with q along two different high- 
symmetry directions, namely T — M [panel (a)] and T — K 
[panel (b)] is shown. It is clear that Im[g(q, cj)] exhibits 
a highly isotropic character. In all the carried tests the 
same isotropic behaviour of the surface loss function was 
found independently of the film thickness. Thus, from 
here on only results for q along T — M are shown, as the 
grid used in this high-symmetry direction is finer than 
the one along T — K. 
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FIG. 2: (Color online) Calculated surface loss function for 
the 3MLs thick Pb(lll) film, with q along two different high- 
symmetry directions, T — M (a) and T — K (b) 
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FIG. 3: (Color online) Surface loss function surfaces Im[g(q, cj)] for 1-8 MLs thick Pb(lll) films calculated at the scalar- 
relativist ic level. The momentum transfer vector q is along T — M. The colour bar applies to all plots, with its maximum 
normalized individually for each case. Solid lines represent the dispersion of the classical modes, see Eq. ([l]). 



B. General results 



The general results of the present work are shown in 
Fig.[3J In order to get insight, the dispersion of the classi- 
cal modes uj± = oj±{q) given by Eq. ^ is represented by 
solid green curves. For each freestanding slab, uj±{q) are 
plotted for an effective thickness corresponding to a num- 
ber of interlayer distances equal to the number of MLs 
forming the slab, as the jellium edge in the first-principles 
calculations was fixed at half an interlayer distance away 
from the outermost atomic layers. 

The results of the present work as plotted in Fig. |3] 
show several modes of different character. First, the 
low-energy symmetric mode is detected for the thinnest 
slabs at small momentum transfer values, closely follow- 
ing the dispersion described by the low-energy U- mode 
of Eq. (IT]) for all thicknesses as represented in Fig. Is] by 
the bottom green line in each panel. However, notice 



that it disappears upon entering the almost dispersion- 
less peak present around cj ~ 7 eV for all thicknesses. 

Also, the high-energy plasmon mode analogous to the 
classical thin film a;+ mode is found for thicknesses 
greater than 2MLs. Note that it is placed at too high en- 
ergies in comparison with the predictions of Eq. (IT]) . Un- 
fortunately, calculations including the 5d semicore elec- 
trons are too computationally demanding in the supercell 
scheme used here. Thus, the influence of the semicore 
electrons on Im[g(q, a;)] has not been checked . The in- 
clusion of the polarizable 5d semicore electrons through 
the use of a model dielectric function ea^^ is ambiguous 
and its use has been discarded in this present study. Note 
the 5d electrons have been shown to play a crucial role in 
the dielectric response of bulk Pb, more precisely in the 
optics and dynamics of the main bulk plasmon.^ 

In Ref. I20] using jellium calculations it was shown that 
the antisymmetric mode disappears for film thicknesses 
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FIG. 4: (Color online) Real part of the induced density at g=0.014 a.u. for the 4 MLs Pb(lll) film, as a function of the energy 
transfer uj and the z coordinate. Panel (b) is a zoom in the {uj^z) range delimited by the orange rectangle in panel (a). The 
orange circle in (b) highlights the fingerprint of the low-energy symmetric mode (see the text). Black solid (red dashed) lines 
mark the position of the atomic layers (jellium edges), uj in eV. 



comparable or smaller than the Fermi wavelength of the 
metal when q ^ {). Instead, peaks corresponding to dis- 
crete interband transitions show up. For Pb, using the 
value rf^ = 2.298, one finds A^^ = 7.52 a.u. which 
is roughly 1.4 times the interlayer distance in Pb(lll) 
films. Thus, the 1 and 2 MLs thick Pb(lll) films present 
electronic effective thicknesses equal to 0.7 and 1.4 times 
A^^, respectively. As seen in Fi g. [3} our results are in 
agreement with the work of Ref. I2Q| as far as the disap- 
pereance of the high-energy mode for thin films is con- 
cerned. In the surface loss function of the single mono- 
layer shown in Fig. |3J a manifold of interband peaks is 
present for energy transfers cj > 11 eV [see also the black 
solid curve in Fig. Is] (a)], where the high-energy mode 
should be present (see the upper green line in the first 
panel of Fig. p|. This is a manifestation of strong QSE 
in the surface- loss function of the single Pb(lll) mono- 
layer. The 2 ML Im[g(q, cj)] results (see Fig. [3| corre- 
spond to the transition between the two different thick- 
ness regimes, at L c^ 1.4A^^. 

An important conclusion of the present work is the 
large difference in spectral weight between the low- and 
high-energy modes of the film, in sharp contrast with 
the results reported in Ref. ^20^ for Ag slabs modeled by 
the jellium approximation. The low-energy mode analo- 
gous to the classical symmetric uo- plasmon appears as 
a faint feature in comparison with the rest of the peaks 
present in Im[g(q, co')]. On the contrary, the high-energy 
mode is the most intense feature in the surface loss func- 
tion of freestanding Pb(lll) films, except for the single 



monolayer. In the latter case, a slightly upwards dispers- 
ing interband peak raises at energies uo o:^ 5.5 — 6 eV. 
It exhibits the highest intensity [see Fig. Is] (a) and (b)] 
together with a vanishing linewidth at momentum trans- 
fer smaller than 0.1 a.u. This long-living mode stems 
from transitions between the highest occupied and low- 
est unoccupied QWSs around the SBZ center (see Fig.fTl), 
representing strong QSE. Once more the 2 MLs results 
represent the crossover with larger thicknesses for which 
the quantization of the states is not reflected in the same 
fashion in the calculated surface loss function. Neverthe- 
less, in the evaluated Im[g(q, cj)] corresponding to the 
2 MLs thick slab, still two interband peaks similar to 
the long-living mode present in the single monolayer are 
found, overlapping with each other. However, their in- 
tensity is greatly decreased in comparison with the cor- 
responding feature of the response of the IML Pb(lll) 
film. 

In Fig. [4] the real part of the two-dimensional 
Fourier transform of the induced density [see Eq. (|3|] 
Re[n^"^(z, q, cj)] for the 4 MLs thick slab is shown. The 
results correspond to a momentum transfer of g = 0.014 
a.u. First, note the antisymmetric distribution of the in- 
duced density with respect to z=0. Second, for spatial 
positions inside the film, several changes of the sign of 
Re[n^"^] are found. In panel (b) results for 4 ^ cj ^ 8 eV 
are zoomed in. Interestingly a sharp change of phase of 
Re[n^"^] can be recognized at a; ::^ 7 eV, signaling about 
the presence of the dispersionless peak seen in Im[g(q, uj)] 
at this energy. This is a general behaviour found at all 



thicknesses at cj :^ 7 eV. 



C. Thickness dependence 



Moreover, the fingerprint of the low-energy mode anal- 
ogous to the classical symmetric uj- plasmon is found, 
as marked by the circle in panel (b) of Fig. Hi As can 
be seen, cj ~ 5 eV is the only energy at which there is 
a noticeable weight of Re[n^"^(z, q, cj)] at the center of 
the slab, 2: = 0. This notable distortion in the general 
antisymmetric distribution of the real part of the two- 
dimensional Fourier transform of the induced density sig- 
nals about the presence of the symmetric plasmon mode. 
This faint but appreciable cj-— like fingerprint has been 
found at all thicknesses in which the symmetric mode 
could be resolved. 



(a) 



o 
o 



CO 
CO 

o 



20- 
10- 

20 



q = (1011 4 a.u. 



r- 'I 



J I V 



(b) 



q = 0.0685 a.u. 




q = 0.1256 a.u. 




q = 0.1828 a.u. 




4 6 8 10 12 14 16 

Energy (eV) 

FIG. 5: (Color online) Surface loss function for different val- 
ues of q along V — M and different film thicknesses. Black 
solid, red dashed, green dashed-dotted and blue dashed- 
dotted-dotted curves represent results for 1, 3, 5, and 8 MLs 
thick Pb(lll) films, respectively. The thick orange (grey) 
solid curve stands for the results deduced from bulk calcula- 
tions without (with) inclusion of the ^d electrons. The vertical 
dashed line marks the classical Pb surface plasmon energy of 
9.57 eV, while the shaded energy interval corresponds to the 
electron energy loss experimental value of 10.6ib0.2 eV.^^ 



In order to relate the results for different thicknesses, 
several cuts of Im[g(q|| , uj)] are plotted in Fig. [5 compar- 
ing the surface loss function of Pb(lll) films of distinct 
thicknesses for the same momentum transfer values of q 
along the T — M high-symmetry direction. 

In panel (a) of Fig. Is] the black curve at a; > 11 eV 
shows the manifold of interband peaks which replaces a 
single high-energy antisymmetric mode for the Pb mono- 
layer. An additional important feature in the surface loss 
function results for the 1 ML slab is the long-living in- 
terband peak found at small momentum transfer, seen at 
UJ = 5.5 eV in panels (a) and (b) of Fig. [5] 

For the 3, 5 and 8 MLs thick films the surface plasmon 
is already present at g=0.1256 a.u., as seen in panel (c) 
of Fig. [5] Note that it presents a remarkably smaller 
intensity than the high-energy antisymmetric mode. On 
the other hand, the low-energy symmetric mode can not 
be seen in the scale of Fig. pj as it is a faint feature (see 
Fig.§. 
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FIG. 6: (Color online) Scalar-relativistic surface loss function 
of the 15 MLs thick Pb(lll) film, q along f - M. The green 
lines stand for the dispersion of the classical thin film modes as 
given by Eq. ([l]) . The pink square represents the experimental 
data of Ref. 55 of c^r^ = 10.6 ±0.2 eV. The same colour code 
as in Fig. [3] is used. 



D. Surface plasmon 

In Fig. [5] the vertical dashed hne marks the classi- 
cal surface plasmon energy ujs = cjp/V^ = yl.br^ , 
which for the averaged valence electron density of bulk 
lead rf^ = 2.298 gives the value uf^ = 9.57 eV. Also, 
the results of the experimental electron energy-loss mea- 
surements of 10.6±0.2 eV^ are represented by the thin 
shaded area. As can be seen, the classical expression 
gives a too low value of the surface plasmon energy by 
about 1 eV. 

On the other hand, in the optical limit {q -^ 0) the 
surface response function can be calculated from the bulk 
dielectric function as^^ ^^ 



g{q^O,uj) = 



s^^^^{q^0^uj)-l 



and thus the surface loss function is 



1' 



Im[g(q -^ O^uj)] (X — Im 



1 



_£buik(q^O,cj) + 1 



(7) 



(8) 



In Fig. [5] the orange (grey) thick solid curve represents 
Im[g(q -^ 0,cj)] calculated using Eq. (Is]) and including 
the 5d electrons in the core (valence). The energy of 
this peak at half width at half maximum (HWHM) is 
of 10.85 (9.3) eV with the semicore electrons excluded 
from (included in) the valence configuration. When us- 
ing the bulk dielectric function obtained without includ- 
ing the 5d electrons, the value retrieved is close to the 
experimental one of 10.6±0.2 eV. However, the agree- 
ment is worsen upon taking the semicore electrons into 
account in the evaluation of s^'^^^{q -^ 0,a;). Note that 
the surface loss function obtained from a bulk calculation 
(without the semicore) through Eq. (|8| is in qualitative 
agreeement with the slab surface plasmon for thicknesses 
greater than 2 MLs at momentum transfer values where 
the modes uj± are uncoupled, see panel (c) in Fig. [sj 
Surprisingly, Im[g((7 -^ O^uj)] calculated from the knowl- 
edge of s^'^^^{q -^ 0, uj) shows a faint peak at 7 eV, mim- 
icking the dispersionless feature which plays the role of 
the short wavelength limit of the symmetric mode uj- 
in the thinnest films. This signals about the bulk-like 
character of the aforementioned dispersionless interband 
mode. 

From Figs. |3] and [6J it seems the surface plasmon dis- 
perses roughly linearly with the momentum transfer. In 
Fig .Jt] the calculated dispersion ujs = 0Js{q), with q along 
r-M, is shown for the 7 ML thick Pb(lll) film. Thicker 
films did not present any remarkable difference in the sur- 
face plasmon dispersion. The values of ujg were evaluated 
at the position of the HWHM. As can be seen, the surface 
plasmon presents a fairly linear dispersion as a function 
of the momentum transfer for g' > 0.1 a.u. 

The straight line in Fig. u^ is the result of fitting uj{q) = 
A-\-B-qfoT q > 0.1 a.u. The obtained values of the fitting 
parameters are A = 10.58 eV and B = 7.35 eV-a.u. It is 
interesting to compare this findings with a simple model 
giving a similar behaviour of ujg = 0Js{q). 
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FIG. 7: (Color online) Surface plasmon dispersion for the 7 
MLs thick Pb(lll) film as a function of q (along T-M). The 
circles represent the calculated values of the HWHM position 
at each q. The red solid line is a linear fit of the computa- 
tional results, while the green dashed one stands for the sur- 
face plasmon dispersion in a hydrodynamic approach of the 
jellium semiinfinite surface (see the text). The shaded orange 
square marks the experimental interval of cj^^^ = 10.6 it 0.2 



In a semiinfinite jellium surface, using the so-called on- 
step hydrodynamic approach,'^^^ the following expres- 
sion for the dispersion of the surface plasmon is found at 
long wavelengths: 






(9) 



where ujp/V2 = vTs^ and P = ^/3/5{vF/2f^^, 
being vp = (97r/4)^/^r~^ the Fermi velocity of a free- 
electron gas of average valence electron density parame- 
ter rs. Using r^^ one gets ujp^/V2 = 9.57 eV and /3^^/2 
= 8.801 eV-a.u. This dispersion is plotted in Fig. [T] as 
a green dashed line, while the orange square shows the 
energy interval for the experimentally determined value 
of uf^^ = 10.6 ± 0.2 eV.^^ The dispersion derived from 
the hydrodynamic approach fails in reproducing a correct 
value for the optical surface plasmon energy (as pointed 
above). N ote th at strictly speaking, Eq. ([9| is valid for 
q <C 2ujp/(3^^^lii the case of lead, this gives the condi- 
tion q <C 1.087 a.u. 

Finally, note it is difficult to deduce a value of 0Js{q -^ 
0) from the present calculations, as the surface plasmon 
disperses with the momentum transfer in contrast to the 
classical picture described by Eq. (fTl). In addition, the 
well-known negative dispersion of the surface plasmon 
as a function of q in the long wavelength limit is not 
retrieved in the present work, as even for the thickest 
film studied (15 MLs) the low- and high- energy modes 
are splitted for the smallest values of q used. 



V. CONCLUSIONS 

In the present work, the surface loss function of thin 
Pb(lll) films has been studied by means of a first- 
principles pseudopotential approach using a supercell 
scheme. 



For 1 and 2 MLs thick films strong QSE have been 
found. The high-energy mode is completely absent in 
the dielectric response of the single monolayer. This is a 
direct consequence of the quantization of the electronic 
states, leading instead to the appearance of discrete in- 
terband transitions in the high-energy range at small mo- 
mentum transfer q [see Fig. |5](a)]. This is in agreement 
with a previous work based on the jellium model. ^^ 

Incorporation of the full 3D ab initio band structure 
also shows a new feature. It does not disperse with the 
momentum transfer for films thicker than 2 MLs, pre- 
senting an energy of uj ^ 7 eV. In practice, this new 
mode plays the role of the classical surface plasmon as 
the long-g' limit of the low-energy thin film mode, as cu- 
disappears upon coupling to the dispersionless peak. To 
the best of our knowledge, this is the first work predicting 
the existence of this new mode as the short wavelength 
limit of the low-energy mode, replacing the role of the 
classical surface plasmon of energy ujs = ujp/\/2. Indeed, 
in Ref.[55la value of 7.2±0.1 eV was reported as the aver- 
age energy of a feature below the surface plasmon energy 
in EELS measurements. We identify this feature as the 
dispersionless mode found in the present ab initio study. 
Surprisingly, the optical surface loss function evaluated 
from bulk calculations [see Eq. ^] also shows a faint 
peak at ~ 7 eV. 



Also, the surface loss function calculated from Eq. ([tI 
is in agreeement with the first-principle results [see FigTJS] 
(c) and (d)]. As regards the surface plasmon dispersion, 
a linear dependence with q has been found in the present 
work. Once its dispersion is fitted to a linear function of 
the momentum transfer, extrapolation of the fitting to 
q ^ gives a value of 10.58 eV, in agreement with the 
experimental u^''^ = 10.6±0.2 eV.^^ 

New electron energy loss spectroscopy measurements 
on Pb(lll) thin films are highly desirable to check the 
present predictions and gain further insight in the dy- 
namics of collective electronic excitations of nanostruc- 
tured systems and the consequences of the quantization 
of the electronic states on them. 
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